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reduction of the annealing time that is proportional to In N, which is much less 
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1. Introduction 

It is generally recognized that quantum annealing (QA) algorithms can find applications 
for solving many practically important problems including optimization of complex 
networks, finding the global minimum of multi-valued functions, and cost minimization. 
Instead of classical annealing algorithm in which the temperature effects are utilized, 
the QA can operate, in principle, at zero temperature. In this case, one first formulates 
the optimization problem in terms of finding the ground state of the effective quantum 
Hermitian Hamiltonian, Ho, which describes a quantum time-evolution of a quantum 
register with N qubits. During the slow and directed time-evolution, the quantum 
system tunnels to the global ground state of this effective Hamiltonian. One of the 
main challenges is to accelerate the speed of QA algorithms, so that the annealing 
time grows not exponentially, but polynomially with the size of the quantum register 

[DEI El HE El HE]. 

Many approaches dealing with finding the ground state of the Hamiltonian, Ho, 
using QA algorithms, can be found in [Tj, El EJ, SI E] (see also references therein). 
To illustrate the idea of QA algorithms, consider the time-dependent Hamiltonian, 
Hit) = Ho + T{t)Hi, where Ho is the Hamiltonian to be optimized and Hi is an 
auxiliary ("initial") Hamiltonian. In order to initiate a non-trivial quantum dynamics 
of a register, the condition, [Ho, Hi] ^ 0, is imposed. The external time-dependent 
field, T(t), decreases from a large enough value to zero during the quantum annealing. 
The ground state of Hi is the initial state, which is assumed to be known. Then, 
the adiabatic theorem guarantees that if T{t) decreases sufficiently slow, the system 
approaches the ground state of Ho, at the end of the quantum annealing. 

One of the requirements of the adiabatic theorem is the existence of a finite gap 
between instantaneous ground state and first excited state of the total Hamiltonian, 
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%(t). However, usually in quantum optimization of iV-qubit models the minimal energy 
gap becomes exponentially small, producing an exponential increase of the annealing 
time [21 HI [9l [Hll [H] . Then, the main problem is on how to speed up the performance 
of the QA algorithms and, at the same time, find the ground state of "Ho- 
rn [12] . we proposed an adiabatic quantum optimization algorithm based on non- 
Hermitian quantum mechanics. Our approach demonstrated that when using the non- 
Hermitian term, Hi, an effective level repulsion takes place for the total Hamiltonian, 
greatly increasing the minimum energy gap. 

In [131 E]> we applied this non-Hermitian quantum annealing (NQA) algorithm 
to Grover's problem of finding a marked item in an unsorted database, and to study 
the transition to the ground state in a 1-dimensional ferromagnetic Ising chain. We 
demonstrated that, for a suitable choice of relaxation and other parameters, the search 
time is proportional to the logarithm of the number of qubits (spins). 

In this paper, we apply our previously developed NQA to study the transitional 
dynamics to the ground state of Ho in a 1-dimensional antiferromagnetic Ising spin chain 
in a time-dependent transverse magnetic field. We also assume that the dissipation 
vanishes at the end of evolution. When the annealing is completed, the system is 
governed by the Hermitian Hamiltonian, Ho- 

The main reason for choosing the antiferromagnetic Ising spin chain is that the 
ground state in this model is much more complicated than in the corresponding 
ferromagnetic chain. Namely, it includes rapid spatial oscillations (correlations) of spins. 
This allows us to demonstrate that our (NQA) algorithm works well for this case, and 
significantly decreases the annealing time. Another reason is to consider the behavior of 
quantum correlation functions (which are more complicated than in the ferromagnetic 
chain) and formation of defects generated in the process of annealing. 

We show that the NQA significantly increases the probability for the system to 
remain in the ground state. In particular, a comparison with the results of the Hermitian 
QA reveals that the NQA reaches the ground state of Ho with much larger probability, 
when we use the same annealing scheme for both. We show that the NQA has a 
complexity of order In N, where N is the number of spins. This is much better than the 
quantum Hermitian adiabatic algorithm which has the complexity for this problem of 
order N 2 . 

This paper is organized as follows. In Sec. II, we introduce a dissipative one- 
dimensional antiferromagnetic Ising system in a transverse magnetic field governed by 
a non-Hermitian Hamiltonian. In Sec. Ill, we study the quench dynamics of our system 
using both analytic and numerical methods. In Sec. IV, we study the correlation 
properties of our system and defects formation at the end of evolution. We conclude 
in Sec. V with a discussion of our results. In the Appendix we present some technical 
details. 
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2. Description of the model 

We consider the 1- dimensional antiferromagnetic Ising chain in a transverse magnetic 
field with dissipation governed by the following non-Hermitian Hamiltonian: 

j N 

H =oJ2 + «+l ~ ma n°n), (1) 

/ n=l 

with the periodic boundary condition, (t n+ i = er 1 . The external magnetic field is 
associated with the parameter, g, the spontaneous decay is described by the parameter, 
5, and = (a* ± icr^)/2 are the spin raising and lowering operators. 

In principle, this model can be realized in a chain of the superconducting phase 
qubits with the tunneling to the continuum mostly from the lowest energy levels. We 
apply the standard Feshbach projection method to obtain the effective non-Hermitian 
Hamiltonian. The motivations for choosing the non-Hermitian relaxation term in Eq. 
(JTJ are discussed in the Conclusion. 

The Hamiltonian in Eq. ([1]) can be diagonalized using the standard Jordan- Wigner 
transformation, following procedures outlined in [HI HSJ [161 El HE] • The Jordan- Wigner 
transformation, which is mapping of a spin- 1/2 system to a system of spinless fermions, 
is given by 

< = 1 - 24c n , (2) 
al = i{cl-c n ) n(l-24c m ), (3) 

a z n = -(c n + cl) J](l-24c m ), (4) 

in which the c n are fermionic operators that satisfy anticommutation relations: 
{clnyCn} = 5 mn and {c m ,c n } = {c^, c^} = 0. Applying these transformations, we 
obtain 

j N 

H = - ( c « c n+i + c{ +1 c ri + c n+1 c n + c\c\ +1 + g(l - 2c\c n ) - 2i8c ] n c n ). (5) 

n=l 

The periodic boundary condition imposed on the spin operators leads to the 
following boundary condition for the fermionic operators: 

c N+1 = -e^'d, (6) 

Afp = J2n=i C n c li being the total number of fermions. From Eq. <^ it follows that 
Mf = N/2 — S x , where S x is the total x-component of the spins. For the particular 
choice of parameters, g = 5 = 0, we obtain, S x = and Mf = N/2. Thus, for S x = 
and Mf — N/2 we obtain periodic (antiperiodic) boundary conditions if N/2 is odd 
(even). Note that since the parity of the fermions is conserved, the imposed boundary 
conditions are valid for all values of the parameters, g and S. 
Next, by applying the Fourier transformations, 

-mt/4 _ 

^^c/^, (7) 

k 
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we find that the Hamiltonian (jSJ) can be recast in Fourier space as follows: 

H = - ^2 2 ( 2 (# ~~ cos ^*)4 c fc _ 9 + shi tp k (c\c ] _ k + c_ fc c fc )) , (8) 

it 

where g = g + i5 and = 2nk/N. For periodic boundary conditions, cjv + i = ci, the 
wave number, k, takes the following discrete values: 
N N 

k = --,...,0,l,...,--l, (9) 

and for antiperiodic boundary conditions, cyv+i = — ci, we obtain 

* = ± i, ± |,...,±^i, (10) 

In what follows, we impose the antiperiodic boundary conditions for the fermionic 
operators. 

The Hamiltonian, H, can be diagonalized using the generalized Bogoliubov 
transformation [18]. Its spectrum is given by £±(k) = Eo ± in which e = 
J cosifik — iJ8, and = J \/ g 2 — 2g cos ^ + 1. There are two (right) eigenstates for 
each k, 

M*)> = ( cos | ), (ii) 

\ sin -j ' 

!«-(*)> = ( " sin ^ | • (12) 



cos y 



in which 



COS(^ fc -^ 

cosgfc = -7== = » (13) 

- 2^( cos + 1 

sm ^ = /~ 2 9 ~ =v (I 4 ) 

- 2 # cos ip k + 1 

with being a complex angle. 

Since for each k, the ground state lies into the two-dimensional Hilbert space 
spanned by |0)fc|0)_fc and |l)fc|l)_fc, it is sufficient to project H k onto this subspace. 
For a given value of k, both of these states can be represented as a point on the complex 
two-dimensional sphere, S%. In this subspace, the Hamiltonian, Hk, takes the form 

H k = e l-j( ~ 9 - C ° S ^ ^ V (15) 
^ sin v?fc ~g + cos (p k J 

Its ground state can be written as a product of qubit-like states, \ip g ) = (£) k \u-(k)), so 
that: 

iVfc) = (g)(cos||0) fc |0)_ fc -sin| |l) fc |l)_ fc ), (16) 

k 

where, |0)fc, is the vacuum state of the mode c^, and \l) k is the excited state: 
|l)* = ct|0> fc . 



6 



For \g\ 3> 1, the ground state is paramagnetic with all spins oriented along the x 
axis. From Eq. f fT3|) we obtain cos 9k — > —1 as \g\ — > oo. From here it follows that, 
\u-(k)) — > ( q j . Thus, the north pole of the complex Bloch sphere corresponds to the 
paramagnetic ground state. On the other hand, when \g\ <C 1 there are two degenerate 
antiferromagnetic ground states with neighboring spins polarized in opposite directions 
along the z-axis. The real part of the complex energy reaches its minimum at the point 
defined by cos 6^ = 1, and, hence, the south pole of the complex sphere is related to the 
pure antiferromagnetic ground state. 



3. Quantum annealing 

In this section, we consider the NQA for the time-dependent Hamiltonian of Eq. (fll) 
written as 



Ur(t)=Ho + Hi{t), (17) 



where 



J N 



2 

n=l 



n u n+l > 



T - 

Ki(t) = g E (9(t)< ~ im)<°«)- (19) 

n=l 

We use with the ground state of the auxiliary Hamiltonian, "Hi(O), as the initial 
state. This state is "paramagnetic" with all spins oriented along the x axis. For g 1 
the Hamiltonian, "H r (0), is dominated by "Hi(O), and the ground state of the total 
Hamiltonian, % T , is determined by the ground state of %i(Q). The term, %i, causes 
quantum tunneling between the eigenstates of the Hamiltonian, At the end of 

the NQA, we obtain, H t (t) = %q. If the quench is sufficiently slow, the adiabatic 
theorem guarantees reaching the ground state of the main Hamiltonian, Ho, at the end 
of computation. 

As shown in Sec. II, the total Hamiltonian, % T {t), in the momentum representation 
splits into a sum of independent terms, % T {t) = J2k^ k ^)- Each %k acts in the two- 
dimensional Hilbert space spanned by \k\) = and \k ) = |0)fc|0)_fc. The 
wavefunction can be written as, \if)(t)) = Ylk li'kit)), where 

\Mt)) = c Q {k,t)\k ) + ci{k,t)\k 1 ). (20) 

Choosing the basis as, k\ = ^ J J and k = ^ ° j , we find that the Hamiltonian, 
Hk(t), projected on this two-dimensional subspace takes the form 

«•(«) = eoMl-jf 9Wr C ° S ** ='7' V (21) 
y smcpk -g(t) + cosifk J 

where £o(t) = Jcosfk — iJS(t) and g(t) = g(t) + i8(t). Further we assume linear 
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dependence of g(t) on time: 

w>-{*7>V* ISr <w 

where 7 = (<? + iS)/r, and g, 5 are real parameters. 

The general wave functions, \ipk) and satisfy the Schrodinger equation and its 
adjoint equation 

i^\tf> k ) =H h (t)\il, k ), (23) 

-£$ k \ = $ fc |ft fc (*). (24) 

Presenting |^fc(t)) linear superposition, 

\M*)) = (Mt)\h) +v k (t)\k 1 ))e i f £0 ® dt , (25) 

and inserting expression fl25|) into Eq. fl23|) . we obtain 

ink = J{{g- cos ip k )u k - simp k v k ), (26) 
iv k = -J(sin tp k u k + (g - cos tp k )v k ). (27) 

The solution can be written in terms of the parabolic cylinder functions, D u (z): 

U k (z k ) = B k D iVk (iz k ) + y/iv k 'A k D^i Vk - 1 (z k ), (28) 
V k (z k ) = A k D^ iVk (z k ) - i^/iv k ~B k D iUk -i(iz k )). (29) 

Here we introduced the new functions: u k {t) = U(z k ), v k {t) = V(z k ) and 

~(j( r - *) - cosv?fc), (30) 




(31) 

27 

The constants, A k and B k , in Eqs. ( 128]) and (|29|) are determined by the initial conditions. 
In the adiabatic basis the wavefunction, \ip k (t)), can be written as follows: 

\Mt)) = (a fc (*)|u_(M)> + A(t)|w+(A:,t)))e^ £0 ^ t . (32) 
From Eqs. ([25]) and (]32]) it follows that 

"few = u k {t) cos— v fc (t)sm— — , (33) 

and 

fa(t) = v k (t) cos ^ + u fc (*) sin ?M. (34) 

We assume that the evolution begins from the ground state. Then, for any k this 
implies, | ck^(0) | = 1 and (3 k (0) = 0. Since cos 9 k (0) = —1, one has, u k (0) = and 
^fc(O) = —1. Using these results, we obtain B k = 0, and the solution of the Schrodinger 
equation can be written as follows: 

U k = A k ^/Ti7 k D^ iUk ^i(z k ), (35) 
V k = A k D_ tUk (z k ), (36) 
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in which A k = (LUjz^O)))- 1 . 

Since for non-Hermitian systems the norm of the wavefunction is not conserved, for 
each k we define the "intrinsic" probability to stay in the ground state as 

P gs (t) - \ ak ^l\ 2 (<rj\ 
Fk {t) ~ \a k (t)\* + \h(t)\ 2 ' { } 

For each k the evolution of the system is independent, and the probability for the whole 
system to stay in the ground state is given by P gs (t) = Yl k Pk S (t). wna t follows, 
we use the asymptotic expansion for Weber functions to calculate the probability, 
P gs = rifc -PjTI 7 ")' to remain in the ground state at the end of evolution, 

For long wavelength modes with \(p k \ <C vr/4, using Eqs. ( 1551) . ( 1561) . we obtain 

w - - ragg • < 38 > 

\y/Th>k~D- iUk -i(z k (T))\ 2 

Using the asymptotic formulas for Weber functions with the large argument values, and 
taking into account that 5 <C g, we obtain, as in fU. 



where 



I _ e -2nMv k 



Kz/ fc = jTn 2 k 2 /(2gN 2 ), (40) 
Sfcftr) = 25Jr/s 2 . (41) 
For the Hermitian QA, 5 = 0, and Eq. (159]) leads to the Landau-Zener formula [T9| [20]. 

P k (r) = l-expf-^P). (42) 



</iV a 

Similar considerations of the short wavelength modes with 3tt/A < \ip k \ < it yields, 
P k s (r) = l + (D(\J\v k \). For short wavelength modes, with 7r/4 < \(p k \ < 7r/2, employing 
the large-order asymptotics for Weber functions, we obtain, P k s (r) = 1 + 0{1/ ^/\v k \). 

Our theoretical predictions are confirmed by numerical calculations performed for 
up to N = 1024 qubits. (See Fig. [TJ) One can observe that while short wavelength 
excitations are essential at the critical point, at the end of evolution their contribution 
to the transition probability from the ground state to the first excited state is negligible. 

Since for each k, the evolution is an independent, the probability of the whole 
system to remain in the ground state at the end of evolution is the product 

^=nw ( 43 ) 

k 

In the long wavelength approximation one can take into account only (p k = n/N, and 
estimate P gs as 

I _ g-27rT/T{) 

~ 1 _ e -27rr/r _|_ e -2jTT / t -Viz 2 (r) ' ^4) 
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Figure 1. (Color online) Left panel. The probability, P® s , to remain in the ground 
state as a function of the scaled time, s — i/r, for the Hermitian QA (8 = 0, J = 0.5, 
g = 10, t = 10 3 , TV = 1024). From bottom to top: k = 1,4,16,32,64. Right panel. 
The probability, P^ s , to remain in the ground state as a function of the scaled time, 
s = t/r, for the NQA (S = 0.25, J = 0.5, g = 10, r = 10 3 , JV = 1024). Blue line 
(k = 1), red line (k = 4), green line (k = 16), orange line (k = 32), black line (k = 64). 



where r = 2gN 2 /(n 2 J) and ^z 2 (r) = 25Jr/g 2 . 

For the Hermitian QA (5 = 0), Eq. (144")) yields the Landau-Zener formula [191 120] 

P, s = l-e- 2 ^ /ro . (45) 

From here it follows that P gs ~ 1, if r > ro. Thus, the computational time for the 
Hermitian QA should be of order iV 2 . 

For the NQA, assuming r ro, we obtain 

P«s ~ ^- 7- (46) 

2vrr 

From here, in the limit of 5 — > 0, we obtain 



P 9 s -> « I- (47) 



2?rr 

The obtained result is expected, as in this case the time of Hermitian annealing, r, is 
small with respect to the characteristic time, r : r r . 
Next, assuming 

_1° »1. 48 
g z 2ttt 

we obtain 

P 9S «l-^-e- 2 ^ 2 . (49) 

As one can see, P gs ~ 1, if the conditions of Eq. PHI) are satisfied. From (14H|) we obtain 
the following rough estimate of the computational time for NQA: r > (g 2 /2JS) \nN. 
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Figure 2. (Color online) Annealing time, r, as function of N. Left panel. Hcrmitian 
QA, S = ( J = 0.5, g = 10). Right panel. Non-Hermitian QA, from top to bottom: 
5 = 0.25,0.5,1. 

In Fig. [2J the annealing time, function of the number of spins is shown for 

the Hermitian QA (left panel) and for the NQA (right panel). The comparison shows 
that for large number of spins (N ~ 1000) and for 5 > 0.25, the annealing time of the 
NQA is ~ 10 3 times smaller than the time for the Hermitian QA. 

These results demonstrate that even moderate dissipation increases the transition 
probability. The characteristic time of non-Hermitian annealing, even for small but 
finite , is determined not only by the number of spins, (as in the Hermitian annealing), 
but mainly by the dissipation rate, . The non-Hermitian quantum annealing has a 
complexity of order , which is significantly better than the corresponding complexity for 
the quantum Hermitian adiabatic algorithm (~ N 2 ). 

4. Correlation functions and defects formation 

During QA the system does not remain in the ground state at all times. At the critical 
point, the quantum system becomes excited, and its final state is determined by the 
number of defects. In the case of the antiferromagnetic Ising chain, the system ends in 
the state such as 

I • • • nmtmtmmtmtmmtmtm, . . ■) (50) 

with neighboring spins polarized on opposite directions along the z-axis and separated 
by walls (defects) in which the polarization of spins has the same orientation. 

The antiferromagnetic spin-spin correlation functions provide information about 
long-range order and defects formation. We restrict ourselves consider of the 
antiferromagnetic correlation function only in the final state, 

1 N 

*(P) = N E«« +P > - «>«+p»> ( 51 ) 

n=l 
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defined at g{r) = 0. 

Using the Jordan- Wigner transformation, the correlation function can be recast as 
follows: 

1 N 

X ^ = jV ^2( B n A n+lB n +2A n+3 . . . Bn+p^An+p), (52) 
n=l 

where A m = 4> + c m and B m = 4, - c m . 

Next, using Wick's theorem, we can express the average in Eq. ( J52l) in terms of 
contractions of pairs: (B m A n ), (B m Bj) and (A m A n ) [21] . When the pairing (B m Bj) = 
and (A m A n ) = for m ^ n, the correlation function can be written as a determinant 
of the Toeplitz matrix, Xrip) = detG, where 

/ Gq GLi . . . G-p+i \ 
G\ Gq ... G_ p+ 2 



G 



(53) 



\ Gp-i Gp-2 • • • Co / 
Its components are defined as follows: 

G p (r) = -Y \^ 2 -\^\ 2 -^t + u l^ e i( m -n )vk _ (54) 



N 

k 



G k = wz I Vl ' ' L 2 V ,.T 9 (55) 



In the thermodynamic limit, N — > oo, we obtain 

1 r (M 2 - \u\ 2 - i(uv* + u*v))e ik *dy 

27T 7_ 7r |^| 2 + |«| 

Computation of the matrix elements of the Toeplitz matrix for the ground state of 
the Hamiltonian, H , yields the following result [2Tj [22j [23j EM] : 

G,-j- 1 ' , * = ° (56) 
fc 1 0, fc = ±l,±2,... V ; 

Using this result, we obtain the correlation function of the antiferromagnetic chain as 

x(p) = (-i) p - 

For the time-dependent problem, the correlation function is not defined by a 
determinant, since in the general case 

{A m A n ) = S mn + 2i^s(3 mn , (57) 

(B m B n ) = -b mn + 2iQ(3 mn , (58) 

where 

Pmn=^- I |7 - (59) 

2m J_ w \u\ z + \v\ z 

In Fig. [3J the results of our numerical calculation for iV = 512 are shown. As can be 
observed (see inset) for the correlation distance, p ^> 1, we can neglect the contribution 
of all ^sj3 mn . This is valid for both Hermitian QA and for NQA, as well. Thus, the 
long-range correlation function can be calculated using the Toeplitz determinant with 

= i r {m) 

2n ]_„ \vr + \ur 
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Figure 3. (Color online) Dependence of 9/3 on p = m — n, 5 = (blue), 5 = 0.25 
(red), 5 = 0.5 (green), 5 = 1 (black) (t = 1000, g = 10, J = 0.5, iV = 512). 



The asymptotics of the Toeplitz determinant can be obtained by applying the Szego 
limit theorem. We use the results of Ref. [25] to obtain for the Hermitian QA (5 = 0) 
the following asymptotic formula for the correlation function: 

Xip) ~ (-l) p exp ( - 0.174|) cos (\J^j + <Po), (61) 

where the phase, v? , depends on the parameter, | = a/ Jr/2g, (for details see Ref. [25]). 

Asymptotically, the correlation function exhibits decaying oscillations. This agrees 
with the results obtained in Ref. [26] for the ferromagnetic Ising chain, where the 
asymptotic behavior of the correlation function is defined by Eq. flBTl) . without the 
factor, (— l) p . For the ferromagnetic Ising chain the oscillatory behavior, corresponds 
to alternating magnetization signs in neighboring ordered domains, separated by kinks 
[251 [26]. For the antiferromagnetic Ising chain the neighboring spins, being polarized 
on opposite directions, are separated by walls (defects) where the polarization of spins 
has the same orientation. In both cases, the Kibble-Zurek (KZ) correlation length 
[26] , £, determines the characteristic domain size as L = ix^\/2ix j ln2. In Fig. HJ 
the results of our numerical simulations are depicted for N = 512 spins and r = 25. 
The theoretical prediction of the domain size, L w 7.5, agrees with the results of the 
numerical calculation. 

The decaying oscillatory behavior of the correlation functions is confirmed by 
our numerical calculations presented in Figs. 0] and [5] One can observe that for 




Figure 4. (Color online) Spin antiferromagnetic correlation function and schematic 
position of spins. Spins, being polarized on opposite directions inside of domain, 
are separated by walls (red arrows) where the polarization of spins has the same 
orientation. The correlation length and domains size are determined by the KZ 
parameter, £. Red diamonds: correlation function calculated for N = 512 spins 
(t = 25, 8 = 0). Solid line: envelope in Eq. (folj) . 




Figure 5. (Color online) Long range behavior of the correlation function at the end 
of annealing. From bottom to top: 5 = 1,0.5,0.25,0 (r = 1000, g = 10, J = 0.5, 
N = 512). Solid lines present the asymptotics given by Eq. ([5T|) . 



NQA the long-range behavior of the correlation functions is much better then for the 
corresponding Hermitian QA, and the characteristic domain size is much longer for NQA 
than L, even for modest values of 5. 
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To evaluate the efficiency of the QA one should calculate the number of defects. The 
computational time of the QA is the time required to achieve the number of defects below 
some acceptable value. The number of defects is equal to the number of quasiparticles 
excited at g = (final state) [16]. It is given by J\f = (iJj T \J\f\il) T ) , where 

1 - 

# = = 2 EC 1 * * 7 ^)- (62) 

k n=l 

In the adiabatic regime we can approximate the average number of defects as, 

Af=l-P k (r), (63) 
with the lowest <pk = tt/N [26J. Using Eq. (jHJ), we obtain 

e ~2TTT/r -nz 2 (r) 

^ ~ I _ e -27rr/r _|_ e -2nr /t q -^z 2 (t) ' ^4) 

When r ^> r we find 

jV^expf-— Y (65) 

V T J 

So, for large annealing times, the NQA does not have any advantage in comparison with 

the Hermitian QA. 

In the opposite limit, r<T , and under the condition 

2J5r r , , 

>! 66 

g z 2ttt 
the average number of defects is approximated by 

A/^-5L e - 2 ^ 2 . (67) 

2lXT 

For the parameters chosen in our numerical simulations: N = 512, J = 0.5, g = 10, and 
for r = 10 3 we obtain J\f 40e~ 105 . Thus, while for the Hermitian QA, the number 
of defects, separating the domains with antiferromagnetic order, is N ~ 40, for the 
NQA with the same protocol, the number of defects decreases as e~ ws . For instance, 
for 5 = 0.25 we obtain J\f m 3. 

Using Eq. (!5li) . we can calculate the density of defects in terms of the correlation 
function as follows: 

n= lim ^ = I(1 + X (1)). (68) 

As shown in the previous section, during the slow evolution only long wavelength modes 
can be excited. So, one can use the Gaussian distribution by replacing simp « and 
cosy?*; ~ 1. In the limit ^2Jt/ g ^> 1, we can employ Eqs. (1391) - f )4T|) to calculate the 
density of defects, 

i r e-^-^dip 
U ~ 7r J 1 - e~ 2n?Ru + e- 2 ^-^ 2 ' 
where iff. — ?■ as N — >■ oo. Performing the integration with 9ftz/ = Jnp 2 /2g and 
3ft^ 2 = 25Jr/g 2 ) we obtain 



n oe - 25Tj/92 $(l-e- 2 ^ 2 ,i,l), (70) 
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Figure 6. (Color online) Density of defects for NQA as a function of the dissipation 
parameter, 5 (t = 10 3 , J = 0.5, g — 10). 



where 




(71) 



denotes the density of defects for the Hermitian LZ problem [16], and $(a;,a,c) is the 
Lerch transcendent [27J. 

In Fig. [6] the density of defects as a function of the decay parameter, S, is shown. 
As one can see, even moderate dissipation greatly decreases the number of defects in 
the system. 

5. Conclusion 

Previously, research on quantum annealing could be conventionally divided by two parts: 
(i) combining classical computers with quantum algorithms and by (ii) building real 
quantum computers. Many schemes and approaches of quantum annealing algorithms 
have been proposed [U [2J, [3j SI O E], [Tl [HI 1211 |29l 130] (see also references therein) . The 
main objective of these publications is to significantly decrease the time of annealing. 
The very popular test models are the one- dimensional Ising spin chains, ferromagnetic 
and antiferromagnetic, which are also useful for practical purposes. In this case, the 
quantum annealing algorithms are used to find the ground state of the chain. 

The approach presented in this paper is related to item (i) above, and was applied 
for the one-dimensional antiferromagnetic spin chain. We have chosen an auxiliary 
Hamiltonian in such a way that the total Hamiltonian is non-Hermitian. This allows us 
to shift the minimal gap in the energy spectrum in the complex plane, and to significantly 
reduce the time required to find the ground state. Our approach leads to an annealing 
time proportional to: IniV, where N is the number of spins. This is much less than 
the time of Hermitian annealing (~ iV 2 ) for the same problem. We also demonstrated 
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the behavior of quantum correlation functions and their relation to the spin defects 
generated during quantum annealing. 

It is interesting to mention here, that similar to the ferromagnetic chain, the main 
contribution to the destruction of the final ground state in the process of quantum 
annealing is produced by the long- wavelength modes, as expected. This is in spite of the 
fact that in the antiferromagnetic chain the ground state represents a rapidly-oscillating 
spatial structure with short wavelengths (opposite to the ferromagnetic chain). 

The dissipative term which we use corresponds to tunneling of the system to its 
own continuum, as usually happens when one applies the Feshbach projection method 
to intrinsic states in nuclear physics and quantum optics. In our case, the intrinsic 
states are the states of the quantum computer (register). So, the probability of our 
quantum computer to survive during the NQA can be small. That is why we use a ratio 
of two probabilities - the ratio of the probability for the system to remain in the ground 
state to the probability of survival of the quantum computer. This relative probability 
(we call it "intrinsic" probability) is well-defined, and remains finite during the NQA. 
So, our approach cannot be directly used in the experiments on QA, but rather as a 
combination of classical computer and protocols for NQA to significantly decrease the 
time of annealing. Also, the dissipative term which we use in this paper is rather artificial 
in the sense that it has no direct relation to physical dissipative mechanisms. At the 
same time, we note that our dissipation term corresponds, in principal, to tunneling 
effects in the superconducting phase qubits if tunneling is realized mostly from the 
lowest energy levels. 

Our hope is that the observed reduction of the annealing time can serve as a guide 
to finding solutions of the NP-complete problems by using a combination of classical 
computers and NQA algorithms. We also would like to mention that many important 
problems still remain to be considered. One of them is the application of this dissipative 
approach to more complicated Ising models with frustrated interactions, in which the 
ground state can be strongly inhomogeneous. Another approach ( "Ising machine" ) was 
proposed recently in [3U [32] in which both energy dissipation and pumping effects were 
utilized for mapping of Ising models onto laser systems. In this case, the system can 
work in the regime of a limit cycle which could be advantageous for performing quantum 
algorithms. Both of these approaches to applying NQA algorithms to Ising models are 
now in progress. 
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Appendix A. 

Exact solution of the Non-Hermitian Landau-Zener problem 

The non-Hermitian Hamiltonian, % k {t), projected on the two-dimensional subspace 
spanned by \ki) = y J J and \k ) = ^ ° J, takes the form 

«»(*) = ^.Wi - j f sw r C0Sf "- "7' ), (a.i) 

\ sm^fe -#(t) + cosy? fc J 

where e (t) = J cosip k + iJ5(t) and g(t) = g(t) — iS(t). We assume a linear dependence 
of the function, g(t), on time: 

where, 7 = (g — iS)/r, and g, 5 are real parameters. 

The general wave functions, \ipk) and (i/j k \, satisfy the Schrodinger equation and its 
adjoint equation 

i^ k ) =n k (t)\ip k ), (a.3) 

-i^$k\ = $k\U k {t). (A.4) 
Presenting \ip k (t)) linear superposition 

\M*)) = (Mt)\ko) + v k (t)\k 1 ))e i J^ d \ (A5) 
and inserting (1A.5I) into Eq. flA.31) . we obtain 

iu k = J((g- cos ip k )u k - sin v k ) , (A. 6) 

n; fc = -J(sin(p k u k + (g - cosy? fe ) v k ). (A. 7) 

Let z k (t) = e m l A \j2J ~ t) ~ cosip k ) be a new variable. Then, for new 
functions, u k (t) = U k (z k ) and v k (t) = V k (z k ), we obtain 

^-U k = ^U k - Vw~ k V k , (A.8) 
dz k 2 

J-Vk = ~^V k - V^~ k U k , (A.9) 
az k 2 

where v k = ( J/27) sin 2 tp k , and the complex 'time' z k runs fromz fc (0) = e W4^2J/7( 7 r- 
cosy?fc) to z fc (r) = -e ln/4 ^2J/'j cos^fc. 

From Eqs. (I A . 8 j) . (I A . 9 j) . we obtain the second order Weber equations 

j£f^(\ + z -i + lVk ) Uk = a ' (A ' 10) 

-**) v > =o - (A - n) 

Solution of the Weber's equations is given by the parabolic cylinder functions, D_ iUk (±2) 
and D il>k _i(±iz). 
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We obtain the solutions of Eqs. (IA.8[ IA.9[) in the form 

U k (z k ) = B k D iUk (iz k ) + y/iui t A k D-. iVk - 1 (zk), 
V k (zk) = A k D^ iVk (z k ) - i^fiv~ k B k D iVk _ 1 (iz k )), 



(A.12) 
(A.13) 



where the constants, A k and B kl are determined from the initial conditions. 
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